Breaking spaces and forms for the DPG method and applications including Maxwell equations

Discontinuous Petrov Galerkin (DPG) methods are made easily implementable using `broken' test spaces, i.e., spaces of functions with no continuity constraints across mesh element interfaces. Broken spaces derivable from a standard exact sequence of first order (unbroken) Sobolev spaces are of particular interest. A characterization of interface spaces that connect the broken spaces to their unbroken counterparts is provided. Stability of certain formulations using the broken spaces can be derived from the stability of analogues that use unbroken spaces. This technique is used to provide a complete error analysis of DPG methods for Maxwell equations with perfect electric boundary conditions. The technique also permits considerable simplifications of previous analyses of DPG methods for other equations. Reliability and efficiency estimates for an error indicator also follow. Finally, the equivalence of stability for various formulations of the same Maxwell problem is proved, including the strong form, the ultraweak form, and a spectrum of forms in between

A Computational Study of the Weak Galerkin Method for Second-Order Elliptic Equations

The weak Galerkin finite element method is a novel numerical method that was first proposed and analyzed by Wang and Ye for general second order elliptic problems on triangular meshes. The goal of this paper is to conduct a computational investigation for the weak Galerkin method for various model problems with more general finite element partitions. The numerical results confirm the theory established by Wang and Ye. The results also indicate that the weak Galerkin method is efficient, robust, and reliable in scientific computing.Comment: 19 page

A note on symmetry boundary conditions in finite element methods

AbstractA short discussion on two kinds of symmetry boundary conditions in the context of variational formulations and finite element methods is presented. Applications including modeling of open boundaries in computational fluid mechanics are discussed and a numerical example is presented

A Primal DPG Method Without a First-Order Reformulation

We show that it is possible to apply the DPG methodology without reformulating a second-order boundary value problem into a first-order system, by considering the simple example of the Poisson equation. The result is a new weak formulation and a new DPG method for the Poisson equation, which has no numerical trace variable, but has a numerical flux approximation on the element interfaces, in addition to the primal interior variable

Structure, shear resistance and interaction with point defects of interfaces in Cu–Nb nanocomposites synthesized by severe plastic deformation

Atomistic modeling is used to investigate the shear resistance and interaction with point defects of a Cu–Nb interface found in nanocomposites synthesized by severe plastic deformation. The shear resistance of this interface is highly anisotropic: in one direction shearing occurs at stresses <1200 MPa, while in the other it does not occur at all. The binding energy of vacancies, interstitials and He impurities to this interface depends sensitively on the binding location, but there is no point defect delocalization, nor does this interface contain any constitutional defects. These behaviors are markedly dissimilar from a different Cu–Nb interface found in magnetron sputtered composites. The dissimilarities may, however, be explained by quantitative differences in the detailed structure of these two interfaces.MISTI-France Seed Fun

The DPG Method for the Convection-Reaction Problem, Revisited

We study both conforming and non-conforming versions of the practical DPG method for the convection-reaction problem. We determine that the most common approach for DPG stability analysis - construction of a local Fortin operator - is infeasible for the convection-reaction problem. We then develop a line of argument based on a direct proof of discrete stability; we find that employing a polynomial enrichment for the test space does not suffice for this purpose, motivating the introduction of a (two-element) subgrid mesh. The argument combines mathematical analysis with numerical experiments

Equivalence between the DPG method and the Exponential Integrators for linear parabolic problems

The Discontinuous Petrov-Galerkin (DPG) method and the exponential integrators are two well established numerical methods for solving Partial Di fferential Equations (PDEs) and sti ff systems of Ordinary Di fferential Equations (ODEs), respectively. In this work, weapply the DPG method in the time variable for linear parabolic problems and we calculate the optimal test functions analytically. We show that the DPG method in time is equivalent to exponential integrators for the trace variables, which are decoupled from the interior variables. In addition, the DPG optimal test functions allow us to compute the approximated solutions in the time element interiors. This DPG method in time allows to construct a posteriori error estimations in order to perform adaptivity. We generalize this novel DPG-based time-marching scheme to general fi rst order linear systems of ODEs. We show the performance of the proposed method for 1D and 2D + time linear parabolic PDEs after discretizing in space by the nite element method

A DPG-based time-marching scheme for linear hyperbolic problems

The Discontinuous Petrov-Galerkin (DPG) method is a widely employed discretization method for Partial Di fferential Equations (PDEs). In a recent work, we applied the DPG method with optimal test functions for the time integration of transient parabolic PDEs. We showed that the resulting DPG-based time-marching scheme is equivalent to exponential integrators for the trace variables. In this work, we extend the aforementioned method to time-dependent hyperbolic PDEs. For that, we reduce the second order system in time to first order and we calculate the optimal testing analytically. We also relate our method with exponential integrators of Gautschi-type. Finally, we validate our method for 1D/2D + time linear wave equation after semidiscretization in space with a standard Bubnov-Galerkin method. The presented DPG-based time integrator provides expressions for the solution in the element interiors in addition to those on the traces. This allows to design di fferent error estimators to perform adaptivity

Combining DPG in space with DPG time-marching scheme for the transient advection–reaction equation

In this article, we present a general methodology to combine the Discontinuous Petrov-Galerkin (DPG) method in space and time in the context of methods of lines for transient advection-reaction problems. We  rst introduce a semidiscretization in space with a DPG method rede ning the ideas of optimal testing and practicality of the method in this context. Then, we apply the recently developed DPG-based time-marching scheme, which is of exponential-type, to the resulting system of Ordinary Differential Equations (ODEs). We also discuss how to e ciently compute the action of the exponential of the matrix coming from the space semidiscretization without assembling the full matrix. Finally, we verify the proposed method for 1D+time advection-reaction problems showing optimal convergence rates for smooth solutions and more stable results for linear conservation laws comparing to the classical exponential integrators

Energy-norm-based and goal-oriented automatic hp adaptivity for electromagnetics: Application to waveguide Discontinuities

The finite-element method (FEM) enables the use of adapted meshes. The simultaneous combination of h (local variations in element size) and p (local variations in the polynomial order of approximation) refinements, i.e., hp-adaptivity, is the most powerful and flexible type of adaptivity. In this paper, two versions of a fully automatic hp-adaptive FEM for electromagnetics are presented. The first version is based on minimizing the energy-norm of the error. The second, namely the goal oriented strategy, is based on minimizing the error of a given (user-prescribed) quantity of interest. The adaptive strategy delivers exponential convergence rates for the error, even in the presence of singularities. The hp adaptivity is presented in the context of 2-D analysis of H -plane rectangular waveguide discontinuities. Stabilized variational formulations and H(curl) FEM discretizations in terms of quadrilaterals of variable order of approximation supporting anisotropy and hanging nodes are used. Comparison of energy-norm and goal-oriented hp-adaptive strategies in the context of waveguiding problems is provided. Specifically, the scattering parameters of the discontinuity are used as goal